Synchronizing Huygens's clocks 
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We introduce an interaction mechanism between oscillators leading to exact anti-phase and in- 
phase synchronization. This mechanism is applied to the coupling between two nonlinear oscillators 
with a limit cycle in phase space, leading to a simple justification of the anti-phase synchronization 
observed in the Huygens's pendulum clocks experiment. If the two coupled nonlinear oscillators reach 
the anti-phase or the in-phase synchronized oscillatory state, the period of oscillation is different 
from the eigen-periods of the uncoupled oscillators. 

PACS numbers: 05.45.Xt, 45.30.+S 



In 1665, Christiaan Huygens reported the observation 
of the synchronization of two pendulum clocks closely 
hanged on the wall of his workshop, [lj, pp. 357-361]. Af- 
ter synchronization, the clocks swung in exactly the same 
frequency and 180° out of phase. Huygens also noted 
that if the two clocks were hanged in such a way that 
the planes of oscillation of the two pendulums were mu- 
tually perpendicular, then synchronization didn't occur. 
Huygens justified the observed synchronization phenom- 
ena by the "sympathy that cannot be caused by anything 
other than the imperceptible stirring of the air due to the 
motion of the pendulum" , [l[ . 

Recently, Bennett et al., [2|, built an experimental 
device consisting of two interacting pendulum clocks 
hanged on a heavy support, and this support was 
mounted on a low- friction wheeled cart. This device 
moves by the action of the tensions due to the swing 
of the two pendulums, and the interaction between the 
two clocks is caused by the mobility of the heavy base 
of the clocks. With this device, the anti-phase synchro- 
nization mode is reached when the difference between the 
natural or eigen-frequencies of the two clocks is less than 
0.0009 Hz. If the difference between these frequencies is 
larger than 0.0045 Hz, the two clocks don't synchronize, 
running "uncoupled" or in a state of beating death, Q . 
This situation is unsatisfactory when compared with the 
observations of Huygens. For example, a difference of 
order of Aui = 0.0009 Hz for the two pendulum eigen- 
frequencies corresponds to a difference in the lengths of 
the pendulum rods of the order of At = y/gfi^Auj/iir, 
which gives, for I = 1 m and g = 9.8 ms~ 2 , At = 4 mm, 
and for I — 0.178 m (the length of the pendulum rods 
used by Huygens, [l[), A£ = 0.02 mm, a precision that 
Huygens certainly couldn't achieve. According to Ben- 
nett et al. @, p. 578], Huygens's results depended on 
both talent and luck. 

Another experimental model mimicking the Huygens's 
clocks system, consists of two pendulums whose suspen- 
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sion rods are connected by a weak string, and one of the 
two pendulums is driven by an external rotor, Q and 0] • 
In this system, the in-phase synchronization is approxi- 
mately achieved with a small phase shift, and the experi- 
mental measurements and the model analysis both agree. 
The numerical results of Fradkov and Andrievsky for this 
device, [J] , show simultaneous and approximate in-phase 
and anti-phase synchronization, tuned by different initial 
conditions. In another experimental device made of two 
rotors controlled by external torques 

(SH]), Andrievsky 
et al, [5], reported approximate anti-phase and in-phase 
synchronization of the two oscillators. In this experi- 
ment, the synchronization parameter is the stiffness of a 
string connecting the two rotors. 

In these experimental systems, there is no clear evi- 
dence of what mechanism is in the origin of the anti- 
phase synchronization, as described by Huygens. In gen- 
eral, it is believed that if the pendulums have slightly dif- 
ferent periods, the two oscillators may not synchronize, 
0, E|. These experimental studies seem to corroborate 
this conclusion. However, as this special type of collec- 
tive rhythmicity occurs in biological systems and several 
other natural phenomena, where individual periods 
are different, it is important to derive and to understand 
the interaction mechanisms leading to exact synchrony. 

In this paper, we introduce an interaction mechanism 
between oscillators leading to exact anti-phase and in- 
phase synchronization. The coupling between the oscil- 
lators is derived by modeling explicitly the physical pro- 
cesses involved in the interaction. The oscillators under 
analysis can be simple harmonic oscillators, pendulums, 
or nonlinear oscillators with a limit cycle in phase space. 

In the Huygens two pendulum clocks system, the pen- 
dulums are hanged in a common support, and the only 
possible interaction between them is due to the tension 
forces generated by the oscillatory motion of the two 
pendulums. These tension forces propagate through the 
common support, that we consider to be elastic. The role 
of the tension forces in the interaction is corroborated by 
the Huygens's finding that when "the clockfaces were fac- 
ing each other" , 0, p. 359], or the planes of oscillation of 
the two pendulums are mutually perpendicular, no syn- 
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FIG. 1: Model to analyze the synchronization of the Huy- 
gens's two-pendulum clocks system. The two pendulums are 
a representation of two nonlinear oscillators. The interaction 
between the pendulums is done by the tension forces at the 
attachment points, and they actuate through an elastic and 
resistive media. Each attachment points is considered to have 
mass M. 



chronization is observed. In fact, the components of the 
tension forces generated by the motion of the pendulums 
are in the plane of motion of the pendulums. 

To model the Huygens's experiment, we consider the 
geometric arrangement of Fig. [TJ where the two pendu- 
lums have masses m\ and m 2 , and lengths l\ and £2, 
respectively. The pendulums are considered connected 
by a massless string with stiffness constant k. The per- 
turbations that propagate along the string are damped, 
and the damping force is proportional to the velocity of 
the attachment points of the string, with damping con- 
stant p. The string and the damping of the attachment 
points simulate the elasticity and the resistivity of the 
common support of the pendulums. 

We also assume that the attachment points of the pen- 
dulums have equal masses M, and their deviations from 
the rest positions are measured by horizontal coordinates 
x\ and X2 , respectively. As we shall see below, the intro- 
duction of the mass constant M is necessary to obtain 
explicitly the equations of motion. 

The system of Fig. [TJ considered without the damp- 
ing forces, is described by the four degrees of freedom 
Lagrangian, 



x\ 



^m 2 (£l8l + x 2 2 
-l M (±i + xl)- 



2£ 1 x 1 8 1 cos 0i ) -\ 
+ 2£ 2 x 2 02 cos 2 ) 
\k(x 2 - xi) 2 



m\g£\ cos X 
+ m 2 g£2 cos 82 

(1) 

where 81 and 82 are the angular coordinates of the two 
pendulums, g is the acceleration due to the gravity force, 
and the last two terms describe the interaction between 
the two pendulums. From {JJ, t ne Lagrange equations of 
motion of the system of FigUJ are, 



m 1 £ 1 8 1 + fi(8 1 , 81) + migsinfli 
m 2 ^ 2 2 + / 2 (02, 82) + m 2 r/sin0 2 
(M + mi)xi + 2px\ + m\£\8\ cos0i 
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—m\X\ cos tii 
— m 2 x 2 cos 82 
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where we have added the dissipative terms implicit in the 
interaction model of Fig. [TJ p is the damping constant 
of the attachment points, and the functions f 1(81,81) 
and f 2(82, 82) describe the escaping mechanism of the 
clocks. The terms in p, fi and / 2 are dissipative terms, 
not contained in the Lagrangian function UJ . 

The system of equations ([2]) implicitly defines a sys- 
tem of ordinary differential equations. If M > 0, l\ > 0, 
mi > 0, £2 > and m 2 > 0, the system of equations ([2j) 
can be solved algebraically in order to the higher order 
derivatives. Solving the system of equations @ in or- 
der to the higher order derivatives, and introducing the 
assumption of small amplitude of oscillations, we obtain, 



mi*i6|i+/i(0i,ffi) + mi$0i 
m 2 ^202 + 12(82, 82) + m 2 g82 
Mxx- h(8 x ,8 l ) + 2px 1 -m x t 
Mx 2 - f2(82,8\) +2p±2-m2i 



= —m\X\ 

= -m 2 X2 

= k(x 2 - xi) 

= -k(x 2 - Xi) 



(3) 

In the particular case of the pendulum clocks, we con- 
sider that its dynamics is well described by a nonlinear 
oscillator with a limit cycle in phase space, Q. To sim- 
plify, we assume that the individual dynamics of each 
oscillator is described by the second order equation, 



m£8 + f(8; A, 0)0 + mgd = 



where, 
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and A and are positive constants. A simple qualitative 
analysis, shows that the second order differential equa- 
tion ([JJ has a unique limit cycle in phase space, [9(. If 
is small, the small amplitude approximation used in the 
derivation of ((3|) is still valid, and the mean radius of the 
limit cycle in phase space is also small. 

Comparing equation (j4|) with the first two equations in 
©, we take as a model for the interaction between the 
Huygens's clocks the system of equations, 



mi40i + /(0i;Ai, 
m 2 402 + /(02;A 2 , 
Mxi-/(0i;Ai,0i 

MCC 2 -/(02;A2,02 



?i)0i + m\g8\ = —m\X\ 

? 2 )02 + "72.902 = -m 2 x 2 

01 + 2pii — rti\g8\ — k(x2 — x\) 

)0 2 + 2p± 2 - 777 2 g02 = ~k(x 2 - Xi) 

which, by (|5|). is a piecewise linear system of equations 
in an eight-dimensional phase space. 

To simplify further, we analyze the solutions of equa- 
tions ([()]), for the particular parameter values, m\ = rri2 = 
m, £1 = £ 2 = £, Ai = A 2 = A and X = 2 = 0, with 
M > 0. 

In a small neighborhood of the origin in phase space, 
we can add the first two equations in ((6]), and also the 
third and forth equations in ([6]), and we obtain, 
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FIG. 2: Numerical solutions of the system of equations (B), 
with damping function ([5}, describing the coupling of two 
identical pendulum clocks with masses m = 1 and lengths 
£ = 1. The parameters of the damping function are A = 0.1 
and 6 — 0.1. The parameter associated with the coupling 
mechanism of the two-pendulum clocks are: M = 0.1, k = 10 
and p = 0.2. The initial conditions are: #i(0) = 0.2, 
2 (O) = 0.3, zi(0) = 0, 3:2(0) = 0, (9i(0) = 0, 2 (0) = 0, 
ii(0) = and ±2(0) = 0. In a) and b), we show the time 
evolution of the angular coordinates of the two pendulum 
clocks, before and after anti-phase synchronization, respec- 
tively. In c) and d) , we show the asymptotic solutions in a re- 
duced phase space of the two pendulum clocks (c) , and of the 
two attachment points (d). The asymptotic solutions of the 
pairs of coordinates (8i(t), 8i(t)) and (62 (t), @2(t)) coincide. 
In this simulation, the period of oscillations is T = 2.461. For 
comparison, in c), we show the limit cycle solution (thin line 
curve) of the reference equation with damping function 



where, 9 = 9\ + 9 2 , and x = x\ + x 2 . Clearly, the solu- 
tions of the system of equations © are related with the 
solutions of 0, provided, |0i(t)| < 9 and \6 2 {t)\ < 9. 

The two nonlinear oscillators described by © exactly 
synchronize in anti-phase if, asymptotically in time, for 
every h > 0, lim tn _ >00 9 1 (t n ) = - lim tn _ ) . 00 6 2 (t n ), where 
t n = to + nh, n = 0, 1, • • • , and to is the initial time. This 
condition implies that, 



lim {O^tn) 



9 2 (t„))= lim 0(t n ) = O 

t n — >00 



(8) 



So, by ijHJ), if the zero solution of the system of equa- 
tions ||7J| is asymptotically stable, and the fixed points of 
the system of equations © is Lyapunov unstable, then 
the two pendulum clocks synchronize in anti-phase, 
These two stability conditions for the zero fixed points 
of equations © and |(7J) are sufficient to ensure exact 
anti-phase synchronization. 

Writing the system of linear equations ([7]) as a first 



FIG. 3: In-phase synchronization of the two pendulum clocks, 
obtained numerically from the system of equations (|6]). The 
parameters are the same as in Fig. [2] except the damping 
parameter p that, in this case, has the value p = 0.02. The 
initial conditions are the same as in Fig. [2] In this simulation, 
the period of oscillations is T = 0.61. In c), the thin line is 
the limit cycle solution of the reference equation (|4j|, with 
damping function (|5j. 



order system of differential equations, we obtain, 
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where, 9 = £, x = v, a = —g(l + m/M)/£, b = 
2\(l/(£m) + l/(£M)), c = 2p/(£M), d = gm/M and 
e = —2X/M. If the eigenvalues of the characteristic poly- 
nomial of the matrix in ([9|) are non positive, then the 
system of two interacting pendulums clocks have anti- 
phase synchronous solutions, provided the fixed points 
of system © are Lyapunov unstable, [9(. 

The condition of non positivity of the eigenvalues of the 
matrix in ([9]) can be derived from the the Routh-Hurwitz 
criterion. It can be shown (|9() that the eigenvalues of 
the characteristic polynomial of the matrix in ([9]) are non 
positive, provided p > po = (X/£)(l + M/m), and pi < 
p < P2, where p\ and p 2 are the roots of the polynomial, 



P(p) 



imiXp 2 - (4mA 2 + 4MA 2 + gm 3 £)p 
+gm 3 X + 2gm 2 MX + gmM 2 X 



(10) 



To analyze numerically the solutions of the system of 
equations (J6]), for the parameters of the oscillator Q, 
we have chosen the parameter values g = 9.8, m = 1, 
£ = 1, A = 0.1, and 9 = 0.1. In this case, the uncoupled 
nonlinear oscillators have the eigen-period T = 2.008. 

In Fig. [21 we show the time evolution of the two pendu- 
lum clocks starting from two different initial amplitudes 
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FIG. 4: Stable limit cycles in a reduced phase space of the 
asymptotic anti-phase and in-phase synchronized states of the 
two pendulum clocks, for p — 0.02. The small limit cycle cor- 
responds to the anti-phase synchronized state, and the larger 
one to the in-phase synchronized state. The periods of oscilla- 
tions are T an n = 2.463, and Ti„ = 0.61. The two stable states 
are reached changing the initial conditions of the pendulum 
clocks. 

with zero velocity. The coupling parameters are k = 10, 
p = 0.2 and M = 0.1. As p > p — 0.11, and the roots of 
the polynomial (TO]) are p x = 0.121 and p 2 = 24.489. The 
condition of instability of the fixed points of the system 
([6]) has been calculated numerically and is, p < 0.393. In 
the numerical simulations of Fig. [21 after a transient time, 
the exact anti-phase synchronization state is reached. 
The period of the two pendulum clocks is T = 2.461, con- 
trasting with the eigen-period T — 2.008 of the uncoupled 
pendulum clocks. The two pendulum clocks synchronize 
in anti-phase, in a phase space orbit different from the 
one obtained if they were uncoupled. Numerically, the 
anti-phase synchronized state is an isolated closed orbit 
(limit cycle) in the eight-dimensional phase space, and 
the periods of the angular coordinates 0t and of the at- 
tachment points Xi are the same. 

Decreasing the damping parameter p, the exact anti- 
phase synchronization regime still persists below the val- 
ues, pi — 0.121 and po = 0.11. 

Decreasing furthermore the damping parameter p, and 
with the same initial conditions of Fig.[2l for p < 0.06, the 
two pendulum clocks synchronize with the same phase 
(in-phase), Fig. [3] 

Changing the initial conditions in the simulations of 
Fig. El from 6»i(0) = 0.2 and 2 (O) = 0.3, to 6 X (0) = 0.2 
and 82 (0) = —0.3, we obtain a new asymptotic solution, 
and the two oscillators still synchronize in anti-phase. 
For p — 0.02, the system of equations j6|) has two stable 
limit cycles in the eight-dimensional phase space, [9j. In 
Fig. 21 we show the stable limit cycles associated with 
the asymptotic anti-phase and the asymptotic in-phase 
synchronized states of the two pendulum clocks. 

Further numerical analysis shows that there exists a 
transition region with p in the interval [0.06,0.07], such 
that, the anti-phase asymptotic regime corresponds to a 
limit cycle in phase space, and the in-phase regime is as- 
sociated with a quasi-periodic orbit, [9j. This suggests 



that the transition between the anti-phase and the in- 
phase synchronized states is due to a non-local bifurca- 
tion, or a sequence of bifurcations. 

In conclusion, we have proposed a model describing 
qualitatively the anti-phase synchronization of clocks as 
observed by Huygens. The model is consistent with the 
physical mechanism associated with the interactions be- 
tween the nonlinear oscillators. The important new issue 
introduced in the model is the possibility of existence of 
small movements of the attachment points of the pendu- 
lum clocks, a situation clearly avoid in the modern exper- 
imental devices. This explains why modern experiments 
have not been able to reproduce the original Huygen's 
results. Anti-phase and in-phase synchrony is obtained 
with periods different from the eigen-periods of the indi- 
vidual oscillators. This shows that the equality between 
the eigen-periods of the individual oscillators is not re- 
quired to obtain the anti-phase synchronization. Drop- 
ping the small amplitude assumption and changing the 
parameters of the individual oscillators, the results pre- 
sented here are still true, [||. 
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